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This  report  consists  of  papers  that  summarize  **or!<  done  on  this  research  project.  The 
£ "najor  results  include:  1)  The  development  and  application  of  the  renormalization  aroup 
method  to  the  calculation  of  fundamental  constants  of  turbulence,  the  construction  of 
turbulence  transport  models,  and  large-eddy  simulations;  2)  The  application  of  RNG  methods 
to  turbulent  heat  transfer  through  the  entire  ranee  of  experimentally  accessible  Peynblds 
numbers;  3)  The  discovery  that  high  Peyr.olds  number  turbulent  flows  tend  to  act  as  if  they 
nad  weak  nonlinearities,  at  least  when  viewed  in  terms  of  suitable  'cuasi-particles'; 

I)  The  further  analysis  of  secondary  instability  mechanisms  in  free  shear  flews,  including 
dhe  role  of  these  instabilities  in  chaotic,  3-D  free  shear  flews;  5)  The  further  develop¬ 
ment  of  numerical  Simula ticns  of  turbulent  spots  in  wall  bounded  shear  flows;  6)  The 
;tudy  of  cellular  automata  for  the  solution  of  fluid  mechanical  problems;  7)  The  clarifi- 
:ation  of  the  relationship  between  the  hyperscale  instability  of  anisotropic  small-scale 
lew  structures  to  long -wave  length  perturbations  and  the  cellular  automaton  description 
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Positive-  and  negative-effective-viscosity  phenomena  in  isotropic  and  anisotropic  Beltrami  flows 

B.  J.  Bayly  and  V.  Yakhot 

Program  in  Applied  and  Computational  Mathematics,  Princeton  University,  Fine  Hall,  Princeton,  New  Jersey  08544 

(Received  23  December  1985) 

A  field-theoretic  approach,  analogous  to  Kraichnan't  direct-interaction  approximation,  to  the  sta¬ 
bility  theory  of  complex  three-dimensional  flows  is  developed.  The  long-wavelength  stability  of  a 
class  of  Beltrami  flows  in  an  unbounded,  viscous  fluid  is  considered.  We  examine  two  flows  in  de¬ 
tail,  to  illustrate  the  effects  of  strong  isotropy  versus  strong  anisotropy  in  the  basic  flow.  The  effect 
of  the  small-scale  flow  on  the  long-wavelength  perturbations  may  be  interpreted  as  an  effective 
viscosity.  Using  diagrammatic  techniques,  we  construct  the  first-order  smoothing  and  direct- 
interaction  approximations  for  the  perturbation  dynamics.  It  is  argued  that  the  effective  viscosity 
for  the  isotropic  flow  is  always  positive,  and  approaches  a  value  independent  of  the  molecular 
viscosity  in  the  high-Reynolds-number  limit;  this  flow  is  thus  stable  to  long-wavelength  distur¬ 
bances.  The  anisotropic  flow  has  negative  effective  viscosity  for  some  orientations  of  the  distur¬ 
bance,  and  is  therefore  unstable,  when  its  Reynolds  number  exceeds  ''Si. 


I.  INTRODUCTION  the  velocity).  This  means  that  the  nonlinear  term  in  the 

Navier-Stokes  equations  reduces  to  the  gradient  of  a  sca- 
One  of  the  most  intriguing  problems  in  the  theory  of  l»r,  which  can  be  absorbed  into  the  pressure  field.  Thus, 

hydrodynamic  turbulence  is  the  formation  of  large-scale  *ny  velocity  field  of  the  form  (1.1)  is  a  nontrivial  exact 

structures  in  a  liquid  performing  random  (turbulent)  steady  solution  of  the  Euler  equations.  Alternatively,  we 

motion  at  small  scales.  It  is  natural  to  explain  the  spon-  can  think  of  the  flow  (1.1)  as  an  exact  steady  flow  of  a 

taneous  formation  of  large-scale  structures  as  a  manifesta-  viscous  fluid  maintained  by  an  externally  imposed  body 

tion  of  a  long-wavelength  instability  of  the  corresponding  force  f  =»(pv/Z.2)V,  where  v  is  the  kinematic  viscosity  of 

small-scale  flow.  In  a  theoretical  approach  to  the  descrip-  the  fluid. 

tion  of  this  phenomenon,  it  is  convenient  to  assume  that  From  the  class  (1.1)  of  Beltrami  flows,  we  shall  select 
the  small-scale  flow  is  driven  by  an  external  force  field  two  for  explicit  calculations.  To  observe  the  effects  of  ex- 


which  varies  with  respect  to  the  space  coordinates.  A  pre¬ 
vious  paper1  treated  the  long-wavelength  stability  of  two- 
dimensional,  periodic  eddy  systems,  driven  by  appropriate 
steady  forces;  it  was  shown  that  the  large-scale  instability 
is  likely  to  occur  only  if  the  small-scale  flow  is  sufficient¬ 
ly  anisotropic. 

In  this  paper  we  shall  consider  small-scale  flows  of  the 
following  form: 

V(x)~  2  ^(Qxfi+'Qxa)*'^171 ,  (1.1) 

Q€  w 

where  If'  is  a  finite  set  of  unit  vectors  which  contains  -Q 
whenever  it  contains  Q.  The  velocity  field  is  guaranteed 
to  be  real  by  requiring  that  the  Fourier  amplitudes  A 
satisfy  Ai  —  Q )  =  /<*( Q )  for  each  QE  W.  The  vector  d  is 
a  unit  vector  orthogonal  to  Q.  For  our  purposes  it  is  un¬ 
necessary  to  specify  just  which  vector  3  is  chosen  out  of 
the  unit  circle  of  candidates;  changing  the  choice  of  3  is 
equivalent  to  changing  the  phase  of  the  corresponding 
amplitude  A,  which  turns  out  to  be  irrelevant  to  the  long- 
wavelength  stability  results.  IrrL  is  the  fundamental 
period  of  each  Fourier  component. 

The  flows  (1.1)  have  the  particular  property  that  the 
vorticity  vector  is  everywhere  a  constant  scalar  multiple 
of  the  velocity  vector  (we  shall  call  this  the  strong  Bel¬ 
trami  property,  as  opposed  to  the  weak  Beltrami  property, 
in  which  the  vorticity  is  a  nonconstant  scalar  multiple  of 


treme  anisotropy  in  the  basic  state,  we  shall  consider  the 
simple  flow  obtained  when  W  contains  only  two  wave 
vectors  ±Q,  which  can  be  chosen  to  be  (±1,0,0).  With 
the  choice  3  =  (0,1,0)  and  real  amplitude  At 
*  the  flow  is 

V~C/rm§(0,cos(x/L),  —  sinU/L))  ,  (1.2) 

which  has  the  form  of  a  monochromatic,  circularly  polar¬ 
ized  wave.  For  the  opposite  extreme  of  highly  isotropic 
flow,  let  us  choose  W  to  consist  of  a  large  number  S  of 

A 

unit  vectors  Q  distributed  isotropically  over  the  unit 
sphere.  The  complex  amplitudes  >4  ( Q )  are  all  given 
modulus  lfrm,/v/2N ;  the  phases  of  A[Q)  are  immaterial, 
and  may  be  chosen  randomly.  In  either  case,  Ums  is  the 
root-mean-square  fluid  velocity,  and  is  held  fixed  as  V  is 
varied  to  obtain  different  degrees  of  isotropy  Beltrami 
flows  of  this  type  have  been  used  extensively  in  simula¬ 
tions  of  diffusion  in  helical  turbulence,*  but  the  present 
work  appears  to  be  the  first  investigation  of  their  dynami¬ 
cal  properties. 

The  possibility  that  a  small-scale  flow  with  nonzero 
average  helicity  might  be  unstable  to  the  growth  of  long- 
wavelength  disturbances  has  been  discussed  by  a  number 
of  authors.  Moffatt'  has  shown  that  a  certain  ver\  s\m- 
metric  Beltrami  flow  with  six  Founer  components,  known 
as  the  Amol'd-Beltrami-Childress  (ABC)  flow  .  is  ins  iscid 
ly  unstable  to  large-scale  disturbances  with  the  same  sign 
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of  helicity  (i.e.,  the  same  “handedness")  as  the  basic  flow, 
and  Moissev  et  al.6  have  shown  that,  in  a  highly  compres¬ 
sible  fluid,  any  small-scale  turbulent  flow  with  nonzero 
helicity  can  support  a  growing  large-scale  secondary  flow. 
The  purpose  of  the  present  paper  is  to  perform  an  investi¬ 
gation  of  the  possibility  of  long-wave  instability  in 
viscous,  incompressible  flow,  which  is  typically  the  situa¬ 
tion  of  greatest  physical  interest.  It  is  our  suspicion  that 
under  these  conditions,  most  nontrivial  Beltrami  flows  are 
stable,  and  as  such  may  occur  and  persist  in  a  large 
variety  of  turbulent  flows.7'® 

In  this  paper,  we  consider  the  behavior  of  long- 
wavelength  disturbances  to  viscous  Beltrami  flows,  which 
we  assume  are  maintained  by  steady  external  body  forces. 
The  perturbation  is  expanded  in  powers  of  the  Reynolds 
number  of  the  basic  flow,  and  we  evaluate  each  term  in 
the  long-wave  limit.  This  expansion  is  essentially  the 
same  as  the  expansion  for  the  long-wave  behavior  of  the 
averaged  impulse-response  tensor  in  a  fully  turbulent 
flow,  for  which  Kraichnan  originally  developed  the  so- 
called  direct-interaction  approximation  (DIA).  The  DIA, 
as  usually  used  in  turbulence  theory,  consists  of  an  equa¬ 
tion  for  the  impulse-response  tensor,  coupled  with  an 
equation  giving  the  evolution  of  the  velocity-velocity 
correlation  function.  In  this  paper,  the  basic  flow  is 
known  exactly,  so  we  do  not  need  the  auxiliary  equation 
for  the  velocity  autocorrelation.  The  resulting  system  is 
considerably  simpler  to  analyze  than  most  applications  of 
the  DIA.  We  construct  both  the  first-order  smoothing 
and  direct-interaction  approximations  and  examine  their 
implications  for  the  two  flows  described  above. 

In  the  long-wave  approximation,  the  growth  or  decay 
rate  of  a  large-scale  perturbation  is  a  quadratic  function 
of  the  wave  number  of  the  perturbation.  For  this  reason, 
the  effect  of  the  small-scale  flow  on  the  evolution  of 
long-wave  disturbances  may  be  interpreted  as  the  action 
of  a  tensor  effective  viscosity.  In  the  case  of  the  simple 
Beltrami  flow  (1.2)  with  two  Fourier  components,  only 
the  first  nontrivial  diagram  survives  in  the  long-wave  lim¬ 
it  for  any  finite  Reynolds  number,  and  the  first-order 
smoothing  approximation  correctly  determines  the  exact 
evolution  of  large-scale  perturbations.  The  effective 
viscosity  is  anisotropic,  and  is  actually  negative  for  some 
orientations  of  the  perturbation.  It  turns  out  that  large- 
scale  perturbations  of  suitable  orientation  are  capable  of 
growing  on  the  flow  (1.2)  when  its  Reynolds  number 
exceeds  v7  2. 

For  almost  isotropic  Beltrami  flows  at  low  Reynolds 
numbers,  first-order  smoothing  theory  indicates  that  the 
viscous  damping  of  long-wave  disturbances  is  enhanced 
by  the  reaction  of  the  basic  flow,  the  effective  viscosity 
being  isotropic  and  positive  in  this  case.  At  finite  Rey¬ 
nolds  number,  however,  all  terms  in  the  expansion  are 
nominally  of  the  same  order  in  the  long  wave  limit,  and 
first-order  smoothing  no  longer  remains  applicable.  The 
solution  of  the  direct-interaction  approximation  for  these 
flows  indicates  that  the  effective  viscosity  increases  with 
the  Reynolds  number,  eventually  approaching  a  value  in¬ 
dependent  of  the  molecular  viscosity  in  the  limit  of  infi¬ 
nite  Reynolds  number.  This  result  supports  the  conjecture 
that  a  large  class  of  Beltrami  flows  with  nontrivial  topolo¬ 


gies  are  stable,  and  might  be  observed  in  experiments  or 
numerical  simulations. 

II.  FORMULATION  OF  THE  LONG  WAVE 
STABILITY  PROBLEM 

The  equations  of  motion  for  a  viscous  fluid  in  the  pres¬ 
ence  of  a  steady,  quasiperiodic  force  field  are  most  con¬ 
veniently  used  in  dimensionless  form.  Using  L  as  the 
length  scale  and  v/L  as  the  velocity  scale,  the  dimension¬ 
less  Navier-Stokes  equations  take  the  form 

vf  +  Rv-Vv=  -Vp+VJv  +  f  , 


where 

f(x)=  2F($)e'$1,  F(Q)  =  /I(Q)(S+jQxn) , 


and  the  |  /4(Q)|,  now  dimensionless,  have  been  scaled 
with  respect  to  the  root-mean-square  fluid  velocity  Umt. 
R  is  the  Reynolds  number  based  on  Umi,  L,  and  v. 
Henceforth,  unless  explicitly  stated,  we  shall  work  ex¬ 
clusively  with  the  nondimensional  quantities  as  defined 
here. 

If  we  define  the  four-dimensional  Fourier  transforms  of 
space-time— dependent  quantities  by 

a(k,tu)s  /  ^-e-'<k  *-"la(x,f), 

(2jt) 

r  ,  (2.3) 

o(x,r)s  J  d k  <Wk  ,-“',fl( k.cu), 

then  we  can  write  the  Navier-Stokes  equations  as 

r,(k,<u)  =  G°(k,0)/,(k)6(w) 

+  G°(k,a>)(~\iR)Pijl(k) 

X  f  dqdaVj(q,a)vi(k-q,co-ai  .  (2.4) 

Here  G°(k  ,<i))  =  (  —  ia>  +  k2)  1  is  the  so-called  zero-order 
propagator,  and 

PytMskj^-^kj/k^+k^j-k.kt/k2)  . 

The  tensor  P^  incorporates  the  pressure  field  and  the  in¬ 
compressibility  constraint,  allowing  us  to  replace  the  two 
equations  (2.1)  by  only  one  (2.4). 

The  Fourier  transform  of  the  force  is  a  sum  of  6  func¬ 
tions: 

/.<k)=  2  F(<Q)6(k-Q>  .  <2.5! 

0e» 

The  Beltrami  properly  of  the  basic  flow  implies  ihai  the 
nonlinear  term 

P,j,(k)  f  dq//q)/|(k-q'  <2  6 

vanishes  identically,  and  the  basic  flow  therefore  takes  the 
form 

vikl  —  G°(k,0)f(k'6(w)  . 


o 
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an  infinitesimal  perturbation  u/(k ,<u)  to  the  flow 

D/(k,w)=G°(k,0)/i(k)8(<u)  +  ti/{k,a)) 


(2.7) 


and  consider  its  evolution  according  to  the  linearized 
equation 

v,  ( k,<u)  =  G°(  k,<w )( —  iR  )Ptbt  ( k ) 

X  /  </qG°(q,0)/*(q)t»;(k-q,6>) .  (2.8) 

Because  of  the  coupling  between  the  perturbation  and  the 
small-scale  basic  flow,  a  pure  long-wave  perturbation  is 
impossible,  but  we  can  formally  separate  the  perturbation 
into  long-and  short-wave  components: 


h  tnro*  uruv  rtwtfTVinr>g»ir»inir^yTTTfT  \rwv*r"  w 

where 


u/(k,<t))sO  if  |  k  |  >e  , 


i>/(k,G>)sO  if  |  k  |  <e  , 


(2.10) 


for  some  long-wave  cutoff  0<e«l.  We  define  the 
operators  (  )  and  I  —  {  >  as  the  projections  onto  the 
long-  and  short-wave  subspaces,  respectively,  of  the  func¬ 
tion  space  we  are  working  in. 

The  long-  and  short-wave  projections  of  the  equation  of 
motion  (2.8)  are  then 


L 


_ i 

u/ ( k,a) ) = G°( k,<u )( —  iR )Pabe(k )  f  dqG°(q,0)(/b(q)v?(k-q,<o))  ,  (2.11a) 

i>f>(k,<u)=G°(k,w)(-«)Pttfc(k)  /  dqGHOJ/^qhi/lk-q.tu) 

+G°(k,<u)(— i/n/Wk) / dqG°(q,0)[/rf(q)i^(k— q,ft>)  — {/rf(q)t>^(k— q,<u)>] .  (2.11b) 

The  short-wave  component  can  be  expressed  in  terms  of  the  long-wave  component  by  formally  iterating  (2.1  lb),  which 
yields  an  expression  in  the  form  of  a  perturbation  series  in  powers  of  the  Reynolds  number  of  the  basic  flow.  When  this 
series  substituted  into  (2.11a),  we  obtain  an  equation  involving  only  the  long-wave  component  of  the  perturbation,  which 
determines  the  evolution  of  the  entire  perturbation.  The  analysis  of  this  equation  occupies  the  rest  of  the  paper. 


i 
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III.  FIRST-ORDER  SMOOTHING 

The  first-order  smoothing  approximation  consists  of  taking  only  one  iteration  of  the  short-wave  equation  (2.11b)  be¬ 
fore  substituting  the  result  back  into  the  long-wave  equation.  This  amounts  to  taking  the  first  term  in  the  Reynolds- 
number  expansion  of  the  evolution  equation,  so  the  results  will  be  strictly  valid  only  for  small  Reynolds  number.  In  this 
regime,  the  molecular  viscosity  dominates  the  inertial  processes,  so  all  disturbances  are  damped.  Nevertheless,  we  can  ob¬ 
serve  whether  the  small-scale  flow  tends  to  enhance  or  diminish  the  dissipation,  and  this  information  will  give  us  some 
indication  of  the  situation  at  higher  Reynolds  numbers. 

The  first  iteration  of  (2.1  lb)  yields 

ve>(k.0)=G°(k,a)(-iX)Pe+(k)  f  dqG^q.OI/^q^fk-q.a) .  (3.1) 

Substituting  into  (2.11a)  then  gives  an  appropriate  equation  of  motion  for  the  long-wave  component  alone: 
i//(k,a>)  =  G°(k,ft>)(  —iR)Paie(k)  f  dq,Go(q,,0)/*(ql)Go(k-q',<u)( -t'R^tk-q1) 

X  f  dq2G°lq2,0)/d(qi)v^{k—q'  — q2) 

=  -RJG°(k,<u)P^(k)  2  G0(k-Q',«)/,(A(k-Ol)n(Q1)/r/QJ)i’.<(k-Ql-QJ.w)  (3.2) 

q'.Q:£  W 

after  performing  the  q  integrals.  But  ^(k—  Q  Q2,g>)  is  nonzero  only  if  Q  J=  —  Q  1  +  0(e),  which  implies 
0  J  =  — Q  '»  because  there  is  only  a  finite  number  of  distinct  wave  vectors  in  W  and  we  are  interested  in  the  limit  e— 0. 
Therefore, 

t;.<(k,6i)=-RJG0(k,«)Pate(k)  2  G°(k-Q)Pt*(k-Q>id(Q)i2(6w-&e<,-/e^£jLv<(k,a>)  ,  (3.3) 

Qew 

where  we  have  used  the  fact  that  a  helical  vector  F(Q)  =  d(Q)[n  +  i'Qx8]  has  the  property  that 

F>(Q)F<(  — Q)=  |  d(Q)  |  *<8* ,  (3  4. 


regardless  of  the  phase  of  d(Q)  or  the  choice  of  8. 

Now,  we  can  simplify  (3.3)  by  taking  the  limit  of  small  |  k  |  and  using  the  incompressibility  conditions 


C.F,(Q)«0,  A,t>/(k,w)*0  . 


(3.5) 
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(3.15) 
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The  small-scale  flow  in  this  case  therefore  enhances  the 
damping  of  all  long-wave  disturbances  irrespective  of 
wave-vector  direction  or  polarization,  as  would  be  antici¬ 
pated  for  such  an  isotropic  small-scale  flow.  In  particu¬ 
lar,  the  isotropic  nature  of  the  small-scale  flow  has  elim¬ 
inated  the  special  directions  in  which  perturbations  were 
amplified  or  only  weakly  damped. 

IV.  HIGHER-ORDER  THEORY: 

DIAGRAM  EXPANSION  AND  DIA 

Because  of  the  complexity  of  the  algebraic  manipula¬ 
tions,  we  shall  use  diagrammatic  notation  for  analyzing 
the  perturbation  equations  (2.11).  The  procedure  is  stan¬ 
dard,  but  we  shall  present  the  derivations  in  detail  for  the 
sake  of  completeness.  Let  the  thin  line  represent  the  zero 
order  propagator  G°,  the  cross  denote  the  force  /,  and  the 
vertex  the  quantity  (-iR)P,  with  indices  to  be  specified 
later.  Also,  let  the  open  bar  represent  the  long-wave  part 
of  the  perturbation,  and  the  slashed  open  bar  the  short¬ 
wave  part.  These  symbols  are  drawn  in  Figs.  1(a)— 1(e), 
respectively. 

Equations  (2.11a)  and  (2.11b)  can  then  be  represented 
by  the  diagrams  in  Figs.  2(a)  and  2(b).  Each  vertex  has  a 
line  exiting  on  the  left,  a  line  entering  from  the  right,  and 
a  line  entering  from  the  upper  right,  whose  other  end  con¬ 
sists  of  a  force  term.  Each  of  these  lines  possesses  a  wave 
vector  and  an  index,  which  must  agree  with  the  wave  vec¬ 
tor  and  index  at  the  force  or  vertex  at  each  of  its  ends.  At 
each  vertex,  the  wave  vector  on  the  left  must  equal  the 
sum  of  the  right  wave  vector  and  the  slanted  force  wave 
vector.  All  repeated  indices  are  summed  over,  and  all 
internal  wave  vectors  are  integrated,  resulting  in  wave- 
vector  sums  of  the  same  form  as  (2.13).  Since  the  vertex 
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FIG.  1.  Constituents  of  diagrams. 

contains  a  factor  of  the  Reynolds  number,  a  diagram  with 
n  vertices  corresponds  to  a  quantity  proportional  to  R  *. 

The  short-wave  perturbation  may  be  expressed  in  terms 
of  the  long-wave  part  by  iterating  Fig.  2(b).  The  first  two 
approximations  are  shown  in  Fig.  3.  Now,  the  term  a  in 
Fig.  3  translates  into  the  algebraic  expression 

- R 2  2  f>a4e(k)G0(Q,)Ft(Q1)G0(k-Q1) 

6'.QJ€W 

xFc*<k-Ql)G0(Q2)Fd(Q2)u,<(k-Q'-QJ,&)) 
and  so  the  only  contribution  to  the  long-wave  projection  b 

^  *  A  « 

of  term  a  comes  from  pairs  Q  ,Q2  that  cancel.  This  can 
be  indicated  diagrammatically  by  replacing  the  force 
branches  by  an  arc  as  in  Fig.  4(a). 

In  general,  the  long-wave  part  of  a  diagram  with  an  odd 
number  of  force  branches  vanishes,  and  the  long-wave 
part  of  a  diagram  with  an  even  number  of  branches  is  the 
sum  of  all  possible  diagrams  obtained  by  connecting  the 
branches  in  pairs.  For  example,  the  long-wave  parts  of 
the  next  nontrivial  diagram  is  shown  in  Fig.  4(b).  Long¬ 
wave  diagrams  that  can  be  disconnected  by  severing  one 
G°  line  are  called  reducible  diagrams  (labeled  r  in  Fig.  4), 
and  those  which  cannot  (labeled  /)  are  irreducible.  So  as 
not  to  count  diagrams  with  internal  wave  vector  k  twice, 
we  require  the  following  rule  for  evaluating  irreducible  di¬ 
agrams:  the  sum  over  all  wave-vector  combinations  must 
exclude  those  for  which  one  (or  more)  of  the  internal  wave 
vectors  is  k.  So,  for  instance,  the  diagram  in  Fig.  4(c) 
with  the  crossed  arcs  has  the  value 


R*  2  G°(k,ft.)/,afcc(k)Ft(Q1)G°(k-Q1,w)Pcd,(k-Q,)fv(Q2)- 

OV-Q1 

X  G°(k— Q  '—Q2,<u )/*,/, (k-Q  '  — QJ)F/(  -  Q  ')G°(k— Q2,tu)P^B1(k— QJ)F/  —  Q2)u^(k,w)  . 


We  can  now  extend  the  iteration  scheme  for  u^lk.w)  to 
arbitrary  order  in  R.  The  next  approximations  are  shown 
in  Fig.  5.  The  two  reducible  diagrams  in  Fig.  5(b)  cancel, 
leaving  an  expression  for  the  fourth  approximation  that 
involves  only  irreducible  diagrams.  This  cancellation  of 
reducible  diagrams  is  quite  general,  and  the  full  expansion 
of  u^fk.w)  is  given  in  Fig.  6.  The  quantity  Im  is  the  sum 
of  all  irreducible  diagrams  with  2m  vertices,  minus  the 
external  legs. 


Having  found  the  expansion  of  v *  in  terms  of  v*  we 
can  now  substitute  it  into  Fig.  2(a)  to  find  the  equation  of 
motion  for  v*  by  itself.  Only  the  diagrams  in  Fig.  6  with 
an  odd  number  of  branches  contribute,  and  all  the  reduci¬ 
ble  diagrams  cancel  in  the  long-wave  projection,  leaving 
the  expression  in  Fig.  7(a).  Figure  7(a)  can  be  rewritten  in 
terms  of  the  renormalized  vertex  Za4(k,cu),  defined  as  the 
sum  of  all  the  legless  irreducible  diagrams.  The  resulting 
equation  [Fig.  7(b)],  whose  algebraic  translation  is 
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FIG.  2.  Equation  (2.1 1)  in  diagram  form. 


FIG  3.  First  two  iterations  of  Fig  2(h) 
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FIG.  4.  Long-wave  components  of  simple  diagrams. 
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d/  ( k,<u )  *  G°(  k,<o  )2«j  ( k,a>  1^ ( k,w ) ,  (4.1) 

can  then  be  solved  provided  we  can  find  some  sensible  ap- 
proximation  for  X^fk.o). 

Formally,  therefore,  we  have  found  the  exact  form  of 
the  Reynolds-number  expansion  of  the  long-wave  equa¬ 
tion  of  motion.  Strictly  speaking,  the  expansion  pro¬ 
cedure  is  only  valid  when  R  «1,  in  which  case  we  can 
get  a  good  approximation  to  the  dynamics  by  taking  only 
the  first  term  in  the  series.  This  term,  which  consists  of 
the  single  second-order  diagram,  gives  the  first-order 
smoothing  approximation  that  was  discussed  in  the  last 
section.  When  the  Reynolds  number  ceases  to  be  small, 
the  series  in  its  raw,  untreated  form  is  not  likely  to  be  use¬ 
ful,  and  we  must  resort  to  more  sophisticated  methods  of 
series  summation. 

Consider  the  subseries,  exhibited  in  Fig.  8(a),  of  the 
series  defining  the  vertex  function  2.  The  sum  of  the 
inner  parts  of  these  diagrams  formally  defines  a  function 
Gat(k,ai)  called  the  full  propagator,  which  we  denote  by  a 
heavy  line.  Every  diagram  in  Fig.  7(b)  belongs  to  exactly 
one  series  analogous  to  Fig.  8(a)  that  can  be  similarly 
summed,  and  so  the  entire  series  for  2  can  be  rewritten  in 
terms  of  G,  as  in  Fig.  8(b).  The  diagrams  remaining  in 
this  consolidated  series  are  those  whose  arcs  cannot  be  di¬ 
vided  into  two  groups  such  that  no  member  of  one  group 
intersects  any  member  of  the  other. 


(g)  = 

•tc. 

FIG.  6.  Full  expansion  of  short-wave  component. 

The  resummation  of  Fig.  7(b)  has  not  simplified  the 
problem  of  calculating  2  at  all;  rather,  we  have  just 
transformed  the  problem  of  finding  2  into  the  problem  of 
finding  G.  But  now  we  can  perform  the  same  resumma¬ 
tion  procedure  on  the  series  for  G  itself,  and  obtain  a  con¬ 
solidated  series  for  G  in  terms  of  itself,  which  can  be  ex¬ 
pressed  in  closed  form  in  terms  of  2  (Fig.  9).  This  equa¬ 
tion,  whose  algebraic  form  is 

Gge(k,tu)  =  G0(k,«a)8at  +  G0(k,<w)2oA(k,ti>)G*c(k,w) , 
can  be  inverted; 

Gai(k,n»)  =  [G0-,(k,«)-2(k,w)]-1^  .  (4.2) 

With  G  given  by  (4.2),  2  is  now  determined  as  the  solu¬ 
tion  of  a  nonlinear  integral  equation,  known  as  the  Dyson 
equation.  Truncating  the  Dyson  equation  and  solving  the 
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FIG.  3.  Third  and  fourth  iterations  of  Fig.  2(b). 
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FIG.  7.  Dispersion  relation  for  long-wave  component 


ivW 


34 


POSITIVE-  AND  NEGATIVE-EFFECTIVE-VISCOSITY  . . . 


387 


4-  — —  4*  — 6— — 


4-  •  •  • 


© 


<0, 


+ 


+  *  «  • 


<» 


FIG.  8.  Resummation  of  subseries. 


resulting  closed-integral  equation  gives  a  new  approxima¬ 
tion  to  the  vertex,  which  may  have  useful  properties  that 
no  finite  number  of  terms  of  the  original  expansion  had. 

The  DIA  for  2  is  given  by  the  lowest  truncation  of  the 
Dyson  equation 

2a/(k,o)=-/lJPa4t(k) 

x  2  -«***&> 

Q€W 

X<7„,(k— $,*>)/>*,(  k—Q)  (4.3) 


with  Gat,  given  by  (4.2).  By  reexpanding  (4.3)  in  powers 
of  R,  it  turns  out  that  the  DIA  simply  consists  of  sum¬ 
ming  up  the  subset  of  the  full  diagram  series  that  consists 
of  all  diagrams  whose  arcs  do  not  intersect  at  all.  Wheth¬ 
er  or  not  the  DIA  is  a  good  approximation  depends  on  the 
relative  importance  of  the  diagrams  with  less  trivial  topo¬ 
logical  structure;  this  issue  will  be  discussed  in  Sec.  VI. 

V.  RESULTS  OF  HIGHER-ORDER  THEORY 

For  the  simple  Beltrami  flow  (1.2),  the  sum  of  all  the 
diagrams  of  a  given  order  n,  which  is  of  order  R"  by  defi¬ 
nition,  is  also  of  order  k*  in  the  long-wave  limit.  In  this 
limit,  therefore,  only  the  single  second-order  diagram  con¬ 
tributes  to  the  dynamics  of  the  perturbation,  and  hence 


=  ~  +  — <=>— 

FIG.  9.  Dyson  equation  for  G.*(k,w). 


the  first-order  smoothing  approximation  is  exact  for 
large-scale  disturbances  to  (1.2).  The  results  of  Sec.  Ill,  in 
particular  the  prediction  that  the  flow  is  unstable  if 
R  >  t/2,  are  consequently  valid  for  all  finite  Reynolds 
numbers. 

One  way  to  demonstrate  this  is  by  an  appeal  to  classical 
(Orr-Sommcrfeld)  stability  theory.  Because  the  velocity 
field  (1.2)  is  independent  of  the  y  and  z  coordinates,  the 
instability  modes  have  the  exact  form 

v  (x,r  )  =  e  (uos.t-os.^os) . 

where 


kos  =  (0,*2,*3),  c=(0,c2,c3). 


and  Uos^os^os  are  functions  of  x  alone.  Eliminating 
u0s  and  ui0s  from  the  linearized  Navier-Stokes  equations 
yields  an  Orr-Sommerfeld  equation  for  uos(xh 


[t/(x)-c] 


d2u  os 
dx 2 


~  kosuos 


d2U(x ) 
dx2 


uOS 


1 

ikosR 


(5.1) 


where  kos=  |kos|,and 


l/U)  =  kos’V(x)//cos  >  c  =  kos-c/kos 


are  the  components  of  V  and  c  parallel  to  kos. 

In  the  long-wave  limit,  with  finite  Reynolds  number, 
the  quantity  kosR  is  small,  and  (5.1)  can  be  solved  to  give 
uos(x)  as  a  power  series  in  kosR.  But  this  is  exactly  the 
same  procedure  that  we  followed  in  Sec.  IV,  although  here 
we  performed  more  preliminary  simplifications  before  ex¬ 
panding.  The  series  obtained  using  either  formalism  must 
agree  at  each  order  in  R\  therefore  the  sum  of  all  the  nth- 
order  diagrams  in  the  diagram  expansion  must  equal  the 
coefficient  of  Rn  in  the  Orr-Sommerfeld  expansion,  which 
is  proportional  to  k  os-  This  result  can  also  be  demon¬ 
strated  directly  by  keeping  track  of  the  orders  of  succes¬ 
sive  approximations  to  v?  while  iterating  (2.1  lb);  we  omit 
this  calculation  because  the  details  are  much  more  intri¬ 
cate  and  no  more  instructive  than  the  Orr-Sommerfeld 
treatment. 

The  reason  for  the  simplification  of  the  series  appears 
to  be  the  fact  that  in  a  sufficiently  regular  flow,  the  long¬ 
wave  perturbation  modes  also  have  regular  spatial 
behavior.  The  Fourier  components  of  the  perturbation  de¬ 
cay  algebraically  with  wave  number  at  any  Reynolds 
number;  indeed,  they  decay  so  rapidly  that  the  component 
of  the  disturbance  at  any  wave  vector  which  is  the  sum  of 
k  and  a  nonzero  combination  of  two  or  more  basic-flow 
wave  vectors  do  not  participate  in  the  dynamics  There  is 
therefore  no  transfer  of  perturbation  energy  to  higher  and 
higher  wave  numbers,  and  it  is  not  surprising  that  the 
long-wave  evolution  depends  strongly  on  the  molecular 
viscosity,  even  in  the  high-Reynolds-number  limit 

For  flows  with  complex  spatial  structure,  the  perturki 
lion  modes  also  have  very  involved  spatial  structure  In 
fact,  the  perturbation  is  likely  to  become  singular  in  the 
limit  of  infinite  Reynolds  number.  There  is  no  finite  set 
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of  terms  from  the  original  expansion  that  satisfactorily 
describes  perturbations  to  almost-isotropic  flows  at 
moderate  or  high  Reynolds  number.  Even  in  the  long¬ 
wave  limit,  all  diagrams  in  the  series  [Fig.  8(b)]  are  nomi¬ 
nally  of  order  k 2,  which  indicates  that  components  of  the 
perturbation  at  all  length  scales  are  important  in  the 
dynamics  of  the  long-wave  part.  It  is  reasonable  to  expect 
that  a  cascade  may  develop  at  high  Reynolds  numbers, 
with  perturbation  energy  from  the  largest  scales  being 
transferred  to  smaller  and  smaller  scales  until  it  can  be  ef¬ 
ficiently  dissipated  by  viscosity.  The  damping  of  long¬ 
wave  disturbances  would  then  become  independent  of  the 
molecular  viscosity  for  large  Reynolds  numbers,  and  be 
determined  by  the  velocity  and  length  scales  of  the  flow 
field  alone. 

Although  the  DIA  does  not  involve  anything  like  the 
complexity  of  the  full  series  for  2  in  Fig.  7,  it  nonetheless 
includes  some  of  the  structure  of  the  perturbation  at  arbi¬ 
trarily  small  scales.  If  there  is  an  energy  cascade  for  the 
perturbation,  we  may  be  able  to  deduce  some  of  its  quali¬ 
tative  characteristics  by  analyzing  its  DIA  equations.  The 
essence  of  the  problem  is  to  solve  (4.3)  for  the  renormal¬ 
ized  vertex  2Lab{k,o).  Since  the  basic  flow  is  isotropic  (to 
an  arbitrarily  good  approximation),  the  tensors  G  and  2 
have  the  forms 

Gab(k,o))=g(k,o)H&eb~kakb/k2)—ih(k,cl>)caby(ky/k) , 

'lab{k,a))=o(k,a)(hab—kakb/k1)+iT{k,a))eebf{ky/k) , 

where  g,  h,  a,  and  r  are  scalar  functions  that  depend  on 
the  wave  vector  k  only  through  its  modulus  k. 

Since  we  are  primarily  interested  in  the  long-time, 
large-scale  behavior  of  the  Green’s  function,  we  can  sim¬ 
plify  the  analysis  of  the  integral  equations  by  setting  the 
frequency  a  to  zero  within  the  integrals;  the  problem  is 
insensitive  to  this  approximation,  and  the  resulting  simpli¬ 
fication  is  enormous.  Then,  by  extracting  the  symmetric 
and  antisymmetric  parts  of  (4.3),  respectively,  we  can 
reduce  the  Dyson  equation  to  a  pair  of  scalar  integral 
equations  over  the  unit  sphere,  which  can  in  turn  be 
transformed  into  integrals  over  the  cosine  of  the  angle  be¬ 
tween  k  and  Q: 

a(k,0)=-~R2  f  dc— — -C  --[a(k,c)g(q,0) 

1  <T 

+(3(,k,c)h{q, 0)] , 

1  ,  ar,  n_ei>  (52) 

r<k,0)=  -  jR2  f  dc - r — [y(*,c)g(9,0) 

1  ? 

-yrj(ktc)h  (q,0)]  , 

where  the  ’’kernels”  of  the  integral  equations  are  given  by 

q=q(k,c)  =  (l+k2-2kc)wi  , 

a(k,c)  =  2kc-4k1c+2k\  0(k,c)=2kc -2k2  , 

y(k,c)  =  4k2c  —4 k1,  rj(k,c)  =  2k2c  —  2k1  . 

The  integral  equations  (5.3)  are  “closed”  by  solving  Eq. 
(4.2)  for  g  and  h  in  terms  of  a  and  r,  which  yields  the  re¬ 
lations 


g(k,w)  =  [k2  —  a(k,a)\/D(k,a)  , 
h  {k,o)=  —Tik,a))/D(k,o) 
with 

Z>(k,Gri  =  [k2  —  <r(k,cu)]2  —  rtfc.a))2  . 

Equations  (5.2)  and  (5.3)  can  be  solved  using  stam 
numerical  techniques.  The  functions  o  and  r  are  appr 
mated  by  specifying  their  values  at  a  suitable  mes) 
points  on  the  k  axis,  and  linearly  interpolating  the  va 
for  all  other  points.  The  integrals  (5.3'  can  then 
evaluated  at  the  mesh  points  using  a  standard  ten-p 
Gaussian  scheme  for  the  c  integrals.  By  performin 
many-dimensional  Newton  iteration,  we  can  rapidly  I 
the  values  of  a  and  r  at  the  mesh  points  which  are  rej 
duced  by  the  integration.  Successively  finer  meshes  w 
used  until  the  values  for  the  eddy  viscosity  (see  beh 
changed  by  1%  or  less  upon  doubling  the  resolution. 

It  is  easily  checked  that  the  asymptotic  behavior 
o(k,0)  is  quadratic  in  k  as  fc— 0,  while  the  asymptr 
behavior  of  iik, 0)  is  cubic  in  k.  So,  although  the  . 
tisymmetric  part  of  2  is  important  for  the  solution  of 
integral  equation,  it  does  not  enter  explicitly  into  \ 
long-wave  equation.  The  behavior  of  long  wave  p 
turbations  is  determined  by  the  quantity 
=  lim*_ol —  o{k,0)/k2R  ],  which  is  a  nontrivial  functi 
of  the  Reynolds  number  R.  Using  Eqs.  (4.1)  and  (5. 
the  long-wave  dispersion  relation  is  given  by 

-iw=— fcJ[l  +  Rr(R)]  ,  (5.4 

or,  in  the  original  dimensional  units, 

-io>=-*2[v+(t/rTnsL)r(R)]  .  (5.4 

Thus,  for  large  Reynolds  number,  a  long- wave  dist  urban' 
damps  out  as  if  acted  on  by  an  effective  eddy  viscosi 
( UmiL)r(R).  A  plot  of  HR)  as  a  function  of  Reynoh 
number  is  shown  in  Fig.  10;  it  should  be  noted  that  as 
increases,  T  tends  rapidly  to  a  finite  asymptotic  valu 
confirming  the  hypothesis  that  at  high  Reynolds  numbe 
the  exact  value  of  the  molecular  viscosity  is  unimportat 
for  the  dynamics.  The  straight  line  through  the  origin  h; 
slope  jf,  corresponding  to  the  first-order  smoothing  a| 
proximation.  The  curve  HR)  is  tangent  to  this  line  at  th 
origin,  confirming  the  result  that  at  low  Reynolds  num 
bers,  the  DIA  reproduces  first-order  smoothing. 

A  couple  of  remarks  may  be  made  about  extensions  o 
the  theory  beyond  the  case  of  long-wave  linear  perturba 
tions.  The  analysis  of  Sec.  II  does  not  actually  depend  oi 
the  precise  character  of  (  )  as  the  projection  onto  func 
tions  with  support  in  a  neighborhood  of  zero  in  wave 
vector  space;  the  formalism  remains  applicable  if  ( 
denotes  the  projection  onto  functions  with  support  in  ; 
neighborhood  of  any  desired  wave  vector.  W'c  can  there 
fore  use  first-order  smoothing  or  the  DIA  to  investigan 
the  behavior  of  a  perturbation  of  arbitrary  wave  vector, 
although  it  becomes  very  difficult  to  evaluate  the  sums  in 
(3.6)  once  k  and  w  become  order  one.  When  k  becomes 
much  larger  than  unity,  corresponding  to  a  very  small 
scale  perturbation,  the  calculations  again  become  simple 
As  expected,  such  high-wave-vector  perturbations  are  rap 
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FIG.  10.  Eddy  viscosity  factor  T  as  a  function  of  Reynolds 
number  R. 


idly  dissipated  by  molecular  viscosity,  before  the  effect  of 
basic  flow  can  be  felt. 

We  can  also  consider  the  effects  of  weak  nonlinearity, 
for  a  perturbation  of  small  but  finite  magnitude.  By  re¬ 
taining  terms  of  second  order  in  the  perturbation  ampli¬ 
tude  in  the  analysis,  a  quadratically  nonlinear  equation  of 
motion  may  be  obtained.  The  coefficient  of  the  quadratic 
term  can  be  expressed  in  terms  of  nontrivially  connected 
diagrams,  and  it  appears  that  even  at  this  order  the  alge¬ 
bra  would  be  too  complicated  to  yield  any  information 
about  the  nonlinear  dynamics.  Whatever  the  form  of  the 
coefficient,  however,  it  must  be  such  as  to  maintain  the 
Galilean  invariance  of  the  whole  problem;  that  is,  if  a  uni¬ 
form  translation  F.filklbtiu)  is  added  to  the  long-wave 
velocity  u^lk.m),  the  frequency  <u  must  change  to 
a— k-V.  This  requirement  alone  is  sufficient  to  deter¬ 
mine  that  the  quadratically  nonlinear  term  has  exactly  the 
same  form  as  the  advective  term  in  the  basic  Navier- 
Stokes  equations;  all  the  complicated  corrections  must 
cancel.  Therefore,  weakly  nonlinear,  large-scale  perturba¬ 
tions  to  almost-isotropic  Beltrami  flow  satisfy 

3,v<+v<-Vv<w— Vp<+vtffVJv<  , 

(5.5) 

v^v+ii/^iriJ!). 

VI.  DISCUSSION 

We  have  considered  the  long-wave  stability  properties 
of  two  viscous  Beltrami  flow  fields,  which  were  chosen  to 
exhibit  the  features  of  extreme  isotropy  versus  anisotropy. 
To  simplify  the  analysis,  we  assumed  that  the  flows  were 
maintained  by  external  body  forces.  Although  such 
forced  flows  rarely  occur  in  the  real  world,  this  hypothesis 
makes  it  possible  to  construct  the  stability  theory  in  an 
uncomplicated  manner.  If  it  were  desired,  our  analysis 
could  then  be  extended,  using  multiple-scale  techniques,  to 
treat  more  realistic  problems,  such  as  the  development  of 
large-scale  structures  in  slowly  decaying  viscous  flows. 

For  the  simple  Beltrami  flow,  the  first-order  smoothing 
approximation  correctly  describes  the  behavior  of  long¬ 
wave  perturbations,  and  instability  is  predicted  if  the  Rey¬ 
nolds  number  of  the  flow  exceeds  v^2.  For  the  class  of 
almost-isotropic  Beltrami  flows,  both  the  first-order 


smoothing  and  direct-interaction  approximations  indicate 
that  the  reaction  of  the  basic  flow  enhances  the  damping 
of  long-wave  disturbances.  A  t  large  Reynolds  numbers, 
the  DlA  predicts  that  the  damping  rate  should  be  essen¬ 
tially  independent  of  the  molecular  viscosity,  which  agrees 
with  an  intuitive  picture  of  the  dynamics.  Our  results  ap¬ 
pear  to  be  in  some  conflict  with  the  general  conclusions  of 
Moffat t5  and  Moissev  et  al.,b  and  it  is  important  to  dis¬ 
cover  whether  this  conflict  is  real  or  whether  it  vanishes 
upon  closer  inspection  of  the  particular  problems 
analyzed.  It  turns  out  that  our  results  are  completely  con¬ 
sistent  with  those  of  Moissev  et  al.b  for  highly  compressi¬ 
ble  flow,  the  essential  difference  being  the  degree  of 
compressibility  involved;  we  illustrate  this  by  repeating 
the  first-order  smoothing  calculation  for  almost-isotropic 
compressible  Beltrami  flow  in  the  Appendix. 

Moffatt  analyzes  the  Amol’d-Beltrami-Childress  (ABC) 
flow,  which  is  a  Beltrami  flow  in  the  class  (1.1)  with  the 
explicit  form 

i7  =Bsin(y)+Ccos(z) , 

F=Csin(z)  +  i4  cos(x)  , 

w  =  A  sin(x)-f  Bcos(y)  , 

where  A,B,C,  are  real  numbers.  A,B,C  are  often  chosen 
to  be  equal  to  unity  in  studies  of  ABC  flows;9  this  choice 
yields  a  flow  which  is  intermediate  in  isotropy  between 
our  simple  flow  and  almost-isotropic  flows.  An  indication 
of  the  inviscid  stability  or  instability  of  the  flow  may  be 
obtained  by  calculating  Amol'd’s  second  variation  of  the 
kinetic  energy  with  respect  to  isocirculational  displace¬ 
ments  of  the  fluid  particles.  This  functional  is  of  indefi¬ 
nite  sign  for  suitably  aligned  long-wave  perturbations  with 
the  same  helicity  as  the  basic  flow,  and  indeed  such  per¬ 
turbations  grow  with  time  in  the  initial  stages  of  their 
development.5  Moffatt’s  analytical  results  are  supported 
by  a  recent  numerical  investigation  of  the  stability  of 
ABC  flow  in- a  viscous  fluid  at  moderate  Reynolds  num¬ 
ber.10 

Despite  its  apparent  isotropy,  the  ABC  flow  has  impor¬ 
tant  properties  that  distinguish  it  from  the  typical 
almost-isotropic  flows  investigated  here.  The  indefinite¬ 
ness  of  the  Amol’d  functional,  for  example,  is  a  property 
of  the  special  orientation  of  the  basic  flow  wave  vectors 
which  ceases  to  hold  when  there  are  more  than  three  in¬ 
dependent  Fourier  components  in  the  flow.  There  are  also 
indications  that  the  form  of  Moffatt’s  equation  of  motion 
for  long-wave  perturbations,  which  involves  the  third- 
order  differential  operator  VxVJ  (hence  the  sensitivity  to 
the  helicity  of  the  disturbance)  gives  way  to  a  reflection- 
invariant  form  involving  only  a  second-order  operator  for 
flows  with  more  than  three  independent  Fourier  com¬ 
ponents.  The  difference  between  the  stability  properties 
of  the  ABC  flow  and  the  almost-isotropic  flows  may  be  a 
consequence  of  the  difference  between  their  topological 
structures.  Recent  investigations  of  the  streamline  struc¬ 
ture  of  the  ABC  flow9  show  that  large  regions  of  space 
are  occupied  by  streamlines  with  highly  regular  structure, 
while  in  Beltrami  flows  with  icosahedral  symmetry,  for 
example,  or  more  complex  cubic  symmetry,  it  appears 
that  almost  all  streamlines  are  ergodic.”  It  seems  likely 
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that  the  ABC  flow  has  more  in  common  with  our  aniso¬ 
tropic  flow  than  our  almost-isotropic  flows,  as  far  as  its 
stability  properties  are  concerned. 

Although  the  DIA  gives  physically  sensible  results  for 
this  and  other  problems,2,1  it  is  still  an  approximation.  In 
fact,  it  is  not  even  a  rational  approximation,  in  the  sense 
that  there  are  no  indications  that  the  omitted  terms  are 
any  smaller  than  those  retained.  The  most  we  can  hope 
when  using  this  technique  is  that  it  incorporate  the  most 
important  physical  processes  sufficiently  completely  to 
give  answers  that  are  qualitatively  correct.  In  this  prob¬ 
lem,  the  most  important  feature  is  probably  the  ability  of 
the  basic  flow  to  mediate  the  transport  of  perturbation  en¬ 
ergy  from  large  scales  to  small  scales.  In  addition,  the 
equations  have  the  property  that  the  nonlinear  terms  con¬ 
serve  the  energy  and  helicity  in  the  flow. 

We  can  actually  construct  a  model  system  with  these 
essential  ingredients.  Imagine  a  large  number  of  realiza¬ 
tions  of  almost-isotropic  flows,  all  chosen  from  an  ensem¬ 
ble  with  well-defined  statistics.  We  can  then  ask  what 
happens  if  each  realization  is  given  a  perturbation  and 
then  allowed  to  evolve  according  to  the  linearized 
Navier-Stokes  equations;  in  particular,  we  want  to  know 
what  happens  to  the  perturbation  averaged  over  all  reali¬ 
zations.  This  situation  is  exactly  analogous  to  that  of  the 
random  harmonic  oscillator  extensively  analyzed  by 
Kraichnan.12 

Now,  consider  the  effect  of  formally  coupling  the 
dynamics  of  the  different  realizations  together.  The 
behavior  will  be  altered,  of  course,  with  the  details  de¬ 
pending  on  the  nature  of  the  coupling.  If  the  couplings 
are  suitably  restricted,  however,  the  total  energy  and  heli¬ 
city  of  all  the  realizations  taken  together  will  still  be  con¬ 
stants  of  motion.  The  system  obtained  by  choosing  the 
couplings  to  be  as  random  as  possible  within  these  restric¬ 
tions  is  called  the  random-coupling  model,  and  its  statisti¬ 
cal  behavior  is  given  exactly  by  the  DIA.12  These  proper¬ 
ties  of  the  random-coupling  model  indicate  that  there  are 
grounds  for  expecting  the  DIA  to  include  at  least  some  of 
the  important  physical  effects  in  other  systems  also, 
indeed,  for  passive-scalar  dispersal  in  reflectionally  invari¬ 
ant  turbulence  the  DIA  and  modifications  thereof  have 
been  highly  successful  in  predicting  the  asymptotic  long¬ 
time  diffusion  rate.2, 4,11 

Unfortunately,  the  DIA  does  not  always  perform  so 
successfully,  particularly  in  the  presence  of  helical  effects. 
The  predictions  of  passive-scalar  dispersal  in  helical  tur¬ 
bulence  are  in  error  by  a  factor  of  three,1  and  helicity  fluc¬ 
tuations  in  conducting  fluids  give  rise  to  anomalous  ef¬ 
fects  in  the  diffusion  of  magnetic  fields  that  are  inaccessi- 
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ble  by  the  DIA.21  There  are  also  indications  that  hel' 
in  the  flow  introduces  effects  that  do  not  appear  at  an 
nite  order  in  the  perturbation  expansion  4  Under  these 
cumstances,  we  cannot  present  our  findings  as  quani 
tive  facts,  but  nonetheless,  we  believe  that  this  work  i' 
important  step  in  understanding  the  dynamical  proper 
of  helical  flows. 
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APPENDIX:  INSTABILITY 
IN  ISOTROPIC  COMPRESSIBLE  FLOW 

Beltrami  flows,  in  particular  the  family  (1.1),  are  al 
exact  solutions  of  the  Euler  equations  in  an  isentropi 
compressible  fluid,  and  can  be  maintained  by  a  suitab 
external  agency  in  a  viscous  compressible  fluid.  In  ord 
to  investigate  the  effects  of  compressibility  on  the  possib 
generation  of  large-scale  motions,  it  is  simplest  at  first  t 
ignore  the  influence  of  pressure  forces  on  the  dynamics  ( 
the  short-wave  part  of  the  perturbation.  This  approxim; 
tion  is  crucial  to  the  analysis  of  Moissev  et  al.,b  wh 
nonetheless  retain  the  pressure  in  the  long-wave  con 
ponent.  For  consistency,  we  shall  neglect  the  perturbatio 
pressure  in  all  its  manifestations  and  obtain  essentially  th 
same  results. 

If  we  neglect  the  perturbation  pressure  field  (and  dro| 
the  incompressibility  condition,  of  course),  the  analysis  o 
long-wave  disturbances  proceeds  exactly  as  in  Sec.  II.  Tb 
only  difference  is  that  the  tensor  through  which  two  velo 
city  components  interact  no  longer  depends  only  on  then 
sum,  but  on  their  individual  wave  numbers 

F^(kl,k2)=kc'6a*-f •  (Al 

This  tensor  is  no  longer  symmetric  in  b.c,  which  is  the 
crucial  feature  allowing  the  long-wave  instability.  Using 
the  first-order  smoothing  approximation  we  obtain  the 
following  equation,  which  is  analogous  to  (3.3),  for  the 
evolution  of  a  long-wave  perturbation: 


u,<(k,a»)=-R2C0(k,«)u#<(k,w)  2  f  M(QMJ^ic(Q.k-Q)G0(k-Q,a»)(5w-CsCrf-»eM1CJ^(-Q.k1]  <A-' 
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For  an  almost-isotropic  flow,  evaluating  the  right-hand 
side  in  the  long-wave  limit  gives 


3,v<  =  -fR2Vxv<+C>(V2)  , 


(A3' 


u«  (k,a>)  =  jR1G0lk,(o)elbeikiv*lk,o>)+Olk2) , 
or,  in  the  space-time  domain, 


which  predicts  instabilities  in  the  form  of  large-scale  heli¬ 
cal  flows  with  helicity  in  the  same  sense  as  the  basic 
flow.5'6 
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